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Abstract 

The problem of estimating a generic phase-shift experienced by a quantum state 
is addressed for a generally degenerate phase shift operator. The optimal positive 
operator-valued measure is derived along with the optimal input state. Two rel- 
■ evant examples are analyzed: i) a multi-mode phase shift operator for multipath 

interferometry; ii) the two mode heterodyne phase detection. 
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i 1 1 Introduction 

C 
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q-i The problem of estimating the phase shift experienced by a radiation beam 

has been the object of hundreds of studies in the last forty years [1]. The 
problem arises because for a single mode of the electromagnetic field there is 
no selfadjoint operator for the phase. This is due to the semiboundedness of 
the number operator [2,3] which is canonically conjugated to the phase as a 
Fourier-transform pair [4]. The most general and, at the same time, concrete 
approach to the problem of the phase measurement is quantum estimation 
theory [5] , a framework that has become popular only in the last ten years in 
the field of quantum information. The most powerful method for deriving the 
optimal phase measurement was given by Holevo [6] in the covariant case. In 
this way the optimal positive operator-valued measure (POM) for phase es- 
timation has been derived for a single-mode field. Regarding the multi-mode 
case, only little theoretical effort has been spent [3] , mostly devoting attention 
to the Lie algebraic structure for two modes [3,7,8]. For two modes one can 
adopt the difference between their photon numbers as the phase shift opera- 
tor, which thus is no longer bounded from below. This opens the route toward 
an exact phase measurement based on a selfadjoint operator [9], with a con- 
crete experimental setup using unconventional heterodyne detection [10,11]. 
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The problem is however complicated by the (infinite) degeneracy of the shift 
operator, and for this reason the optimal states for this case have never been 
derived. 

In this paper the general problem of estimating the phase shift is addressed 
for any degenerate shift operator with discrete spectrum, either S — Z (un- 
bounded), or S = N (bounded from below), or S = 7L q (bounded), generalizing 
the Holevo method for the covariant estimation problem. We find the optimal 
POM for estimating the phase shift of a state |-0o), and then we optimize 
the state itself. The degeneracy of the shifting operator is removed through a 
simple projection technique. The case of mixed input state, which is generally 
very difficult, is considered in some special situations. Two sections are de- 
voted to the analysis of two relevant examples: one concerning a multi-mode 
phase estimation problem that arises in multi-path interferometry; the other 
involving a shift operator that is the difference between the number of photons 
of two modes, corresponding to unconventional heterodyne detection of the 
phase. 



2 Optimal POM for the phase-shift estimation 



We address the problem of estimating the phase-shift pertaining to the 
unitary transformation 



p+ = e-« H p e i * H (1) 

where H is a self-adjoint operator degenerate on the Hilbert space Ti with 
discrete (un)bounded spectrum S = Z, or S = N, or S = Z q , q > 0, and 
p is a generic initial state (actually in the following we will mostly restrict 
to the pure state case). The estimation problem is posed in the most general 
framework of quantum estimation theory [5] on the basis of a cost function 
C(0*, 0) which weights the errors for the estimate 0* given the true value 0. For 
a given a priori probability density po(<f>) for the true value the estimation 
problem consists in minimizing the average cost 



2tt 2tt 

c = J d<t> Po {<t>) J d<t>*c{<t>*,<t>) P {<t>M) , (2) 



where p(<fr*\4>) is the conditional probability of estimating 0* given the true 
value 0. The average cost is minimized by optimizing the positive operator- 
valued measure (POM) [5] d//(0*) which gives the conditional probability by 
the Born rule 
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(3) 



We consider the general situation in which is a priori uniformly distributed, 
i.e. with probability density po((f>) = l/2-n. Moreover, we want to weight errors 
independently on the value of the phase, but only versus the size of the 
error 0* — 0, so that the cost function becomes an even function of only one 
variable, i.e. C(0*,0) = (7(0* — 0). It follows that also the optimal conditional 
probability will depend only on 0* — 0, and the optimal POM can be obtained 
restricting attention only to phase-covariant POMs, i.e. of the form 



d^) = e-^^^, (4) 

where £ is a positive operator, satisfying the completeness constraints needed 
for the normalization of the POM / o 2?r ci/x(0) = I. In fact, using Eq. (3) and 
the invariance of trace under cyclic permutations one can easily recognize that 
p(0*|0) = p(0* — 0) if and only if <ip(0*) is covariant. Hence the optimization 
problem resorts to finding the best positive operator £ for a given cost function 
C(0) and a generic given state po- As we will see, the POM obtained in this 
way is optimal for a whole class of cost functions and initial states po- Once 
the best POM is obtained, one further optimizes the state p . This resorts to 
solving a linear eigenvalue problem. In fact, the average cost can be written 
as the expectation value of the cost operator C, i.e. 



C = Tt[C Po ] (5) 
where 

C = J rfp(0)C(0) . (6) 

Using the Lagrange multipliers method to account for normalization and mean 
energy one has to minimize the function 

£[ Po ] = Tr[C>o]-ATr[p ] (7) 

which for a pure state l^oX^ol is a quadratic form whose minimum is given 
by the eigenvalue equation 



C\ik) = A|Vo> (8) 

with the Lagrange parameter A playing the role of an eigenvalue. The linear 
problem can be easily extended to account also for finite mean energy. 



3 



In summary, our problem is to minimize the cost C for a given cost function 
C(4>) in Eq. (2). This is done in two steps: i) by optimizing the positive operator 
£ for given generic fixed state p : this will give a POM which is optimal for 
an equivalence class of states S(p ); ii) by further optimizing the state in the 
equivalence class £(po)- Since the original state was arbitrarily chosen, this will 
give the absolute minimum cost and the corresponding set of optimal states 
and POM's. 

The solution of the optimization problem is conveniently posed in the repre- 
sentation where H is diagonal. The operator H is generally degenerate, and 
we will denote by \n) v a choice of (normalized) eigenvectors corresponding 
to eigenvalue n, v being a degeneracy index, and by U n the projector onto 
the corresponding degenerate eigenspace. The problem for an input gener- 
ally mixed state p is too difficult to address: therefore, we focus our at- 
tention on the case of pure state po = \i J o){i J o\, and we will leave some 
general assertions on the mixed state case for the following. The problem 
is restricted to the Hilbert space Tin spanned by the (normalized) vectors 
\n) oc n n |^ ) 7^ with the choice of the arbitrary phases such that (n\ip ) > 0. 
Hence the POM can be chosen of the block diagonal form on H = H\\ <8> H±, 
i.e. dp(4>) = dp\\(4>) © dp±((j)) with dp±(if) any arbitrary POM on H±. For 
the optimization of the POM we consider II n |^o) ^ Vn 6 S, as it is clear 
that the resulting POM will be optimal also for states having zero projec- 
tion for some n e S. In this fashion the problem is reduced to the "canoni- 
cal" phase estimation problem restricted to TC\\: \tpo) — > exp(iH\\<f>)\iJ; ) where 
H\\ = J2nes n \n)(n\ and |t/> ) — J2 n es w n\ n )- Now the problem is to find the 
positive operator £y that minimizes the cost C in Eq. (2). On the \n) basis the 
operator £y is written as 



^11 = H \n)(m\^ nm . (9) 

n,mgS 

For a generic even 27r-periodic function C{4>) = — X^o c ^ cos ^ ^ ne average 
cost is given by 



■y oo 

C = -c --^2ci (4'o\n)(m\4) )£ nm . (10) 

1=1 \n—m\=l 

Positivity of £ implies the generalized Schwartz inequalities 



|£nm| — \J ^nn^mm 1 j (H) 

where the last equality comes from the POM completeness / dp\\{(f)) = ly. One 
can write 
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sign(q) (M n )( m \^o)in m < \(^o\n)\\(m\^ )\ , (12) 

\n— m\=l \n— m\=l 

and the equality is obtained only for £ nm = sign(c|„_ m |) (notice that we chose 
(0o | n) > Vn G S). The minimum cost is 



1 oo 

c = -c --YWi\ E |<^b|n>||H^>l (is) 

1=1 \n— m\=l 

where we put sign(O) = 1, since the cost C is independent of £ nm for c\ n - m \ = 0. 
Notice that positivity of £y is not generally guaranteed for any set of sign(q). 
However, one can easily check that £y > if sign(c| n _ m |) = exp[i7r(e n — e m )], 
e n being any integer valued function of n. In fact, this choice corresponds 
to a unitary transformation of the operator £y optimized with all q > 
VZ > 1 (the parameter Co is irrelevant). The particular choice q > VZ > 1 
has been considered by Holevo [6], and includes a large class of cost func- 
tions corresponding to the most popular optimization criteria, as: i) the likeli- 
hood criterion for C(<f>) = — ^(0); H) the 27r-periodic "variance" for C(0) = 
4sin 2 (0/2); Hi) the fidelity optimization C(<f>) = 1 - K^'ole^lV'o)! 2 (here 
q = 2J2\n-m\=i \ w n\ 2 \u! m \ 2 ). For the Holevo class of cost functions the optimal 
POM becomes 



^||(0) = ^|e(0))(e(0)|, (14) 
where the (Dirac) normalizable vectors \e(<f))) are given by 

|e(0)) = J2 e in< "\n) ■ (15) 

The vectors |e(0)) generalize the Susskind-Glogower representation \e 1 ^) = 
Y^=o ein ^\ n ) f° r generic integer spectrum. Therefore, the optimal POM dfi(<f)) 
is the projector on the state \e(4>)) in the Hilbert space Ti\\, and it is orthogonal 
for either S = Z, or S = Z q , whereas it is not for S — N. Notice that the POM 
(14) is also optimal for a density matrix po which is a mixture of states in 7i\\, 
with the additional constraint of having constant phase along the diagonals. 
This can be easily proved by re-phasing the basis \n) in such a way that all 
matrix elements of po become positive. Then the assertion easily follows in a 
way similar to the derivation from Eq. (10) to Eq. (13). Moreover, it is easy to 
see that the pure state case minimizes the cost, which for the optimal POM is 
given by C = — q J2n£s( n \Po\ n + (remember that po > implies that 
|(n|po| m )| 2 ^ (^IpoI^X^IpoI 77 " 1 ), an d the bound is achieved by the pure state 
case (n\po\m) = w^w m ). Finally we want to emphasize that for the bounded 
spectrum S = Z q there is no need for considering a continuous phase dp((p). In 
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fact, it is easy to show [12] that the same average cost is achieved by restricting 
to the set of discrete values {<p s = — , s G Z 9 }, (q = dim(7Y||)), and using 
as the optimal POM the orthogonal projector- valued operator \e(4> s )){e(4> s )\. 

Once the form of the optimal POM is fixed, one can optimize the state |t/> ) 
solving the linear problem in Eq. (8). In the following we show two exam- 
ples of estimation of the phase shift pertaining to highly degenerate integer 
operators (finite dimensional cases are considered in Ref. [12]). In the first 
example we consider the operator H = Y^iLi I (4 a i that describes a multipath 
interferometer, involving M different modes of radiation. In the second, we 
focus our attention on the two-mode phase estimation using unconventional 
heterodyne detection, where the phase shift operator H = a^a — b^b is given 
by the difference of photon numbers of the two modes. 



3 Optimal POM for multipath interferometer 



We consider the operator 



M 



H = ^2 1 a\ai 



(16) 



i=i 



as the generator of the phase shift in Eq. (1). Such phase shift affects a 
M— mode state of radiation in a multipath interferometer, where contiguous 
paths suffer a fixed relative phase shift [13] (this is also a schematic represen- 
tation of the phase shift accumulated by successive reflections in a Fabry-Perot 
cavity). The operator H in Eq. (1) has integer degenerate spectrum S — N. 
We can take into account the degeneracy by renaming the number of photons 
of different modes as follows 



H\n) v = n\n) v , 
with v = (z/ 2 , v 3 , . . . , v M ), and 



(17) 



\n) u = 



M \ 

1=2 I 



\vm) ■ 



The allowed values of v are restricted to the set 8k given by 



V2 = 0, 1, ... , 



A:" 
2 



^3 = 0, 1, 



k — 2z/ 2 
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E^ 1 In 

M 



(19) 



where [x] denotes the integer part of x. 



For the unshifted initial state j^o) we choose a linear symmetrized superposi- 
tion of eigenvectors in Eq. (17), namely 



IV'o) = E w n \n) sym , (20) 



n=0 

where 



|n) sym = -= Y 8 I E^ ~ n ) K) ® W2) <8> ■ ■ ■ ® \vm) , (21) 

V {ui} \i =1 ) 

N n being the number of elements v G S n . Without loss of generality, the basis 
|^)sym has been chosen such that the coefficients w n in Eq. (20) are real and 
positive. According to Eqs. (14) and (15) the optimal POM readily writes as 
follows 



<W) = tt E e * (n ~ m)0 l n )sym sym (m| . (22) 

Z7T n,m=0 

One can now choose a cost function and then minimize the average cost for 
the POM (22) upon varying the coefficients w n of the state (20). By choosing 
the cost function C(<f>) = 4sin 2 (0/2) and by imposing the normalization con- 
straint through the Lagrange multiplier A, the eigenvalue equation (8) gives 
the recursion for the coefficients w n of the form 



w n + w n+2 - 2Aw„+i = . (23) 

The solutions of Eq. (23) can be found in terms of the Chebyshev's polyno- 
mials, and the corresponding optimal state writes as follows 



/2\ 1 /2 00 

W = (-) J2 sin i(n + l)0]\n) sym , = arccos A. (24) 

The state in Eq. (24) is Dirac-normalizable. It is formally equivalent to the 
eigenstate of the cosine operator C of the phase of a single mode [14]. The 
Dirac normalizability comes from the non existence of normalizable states that 
minimize the uncertainty relation for cosine and sine operators 
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AC AS" > - 
~ 2 



<[C,S]> ^^(lOXOl) , (25) 



as proved in Ref. [15]. 



4 Phase-difference of two-mode fields 



In the previous example H was bounded from below and S = N, such that the 
degenerate case is reduced to the standard Holevo's problem. For the difference 
operator H = a^a — b^b one has S = Z, and the set of eigenvectors \d) u can 
be written in terms of the joint eigenvector \n)\m) for the number operators 
a) a and tfb with eigenvalues n and m as follows 

\d) v = \d + v)\v) , 

deZ; v e [max(0, -d),+oo) . (26) 
We consider an initial state \ipo) of the form 



l^o) = ^o|0)|0) + (h n \n)\0) + h_ n \0)\n)) , (27) 

n=l 

where the basis has been chosen to have h n > 0, Vn. The optimal POM writes 
in the form of Eq. (14) in terms of the vectors |A n ), n G Z, where 



IX \ _ J l n )o = \n)\0) , n > , 

|AB '-\|n) w = |0>||n|), n<0. ^ 
Here, the generalized Susskind-Glogower vector \e((f>)) is given by 



+00 

E 

nGZ d=l 



e(0)) = E e m0 |A n ) |0)|0) + £ [f+ |d>|0> + e-^|0)|rf)) . (29) 



Notice that, differently from the usual case of spectrum S — N, now the POM 
is orthogonal (in the Dirac sense): 



+00 



<e(0)|e(0')> = E e**+-n = fatt - 0) , (30) 

n=— 00 

where #271- (0) is the Dirac comb. This means that in this case it is possible to 
define a selfadjoint phase operator 
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+7T 

0=| #|e(0))(e(0)|0, (31) 

— 7T 

as already noticed by Hradil and Shapiro [9,10]. 

We now address the problem of finding the normalized state of the form (27) 
with finite mean photon number that minimizes the average cost evaluated 
through the ideal POM (14). As a cost function we choose again C((f>) = 
4sin 2 (0/2) (periodicized-variance criterion), corresponding to the cost opera- 
tor 



C = 2 - e+ - e~ , (32) 
where 

e+ = £|A n+1 >(A n |, e- = (e+)t. (33) 

Introducing the energy operator E = a) a + tfb and an additional Lagrange 
parameter accounting for finite mean energy (E), the eigenvalue problem in 
Eq. (8) rewrites as follows 



[C-X'-fi'iata + tfb)]^) = 0, (34) 

where A' and ji' are the Lagrange multipliers for normalization and mean 
energy, respectively. The following recursion relations for the coefficients h n is 
obtained 



h n +i + hn-! - n(X+ \n\)h n = , (35) 

with A = (A' — 2)//i' and /i = —/i f . The solution of Eq. (35) is given in terms 
of Bessel functions of the first kind in the following form 



h n = fc(A, /x) J A +|n|(2/» , (36) 
k(X, /i) being the constant of normalization 



+ 0O 



E ^ 2 +H (2//i) 



-1/2 



(37) 



The matching of the recursion for positive and negative indices leads to the 
condition 
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AJ A (2/^) - (2///) J A+ i(2//x) = (2///) 



d 



Ja(2//x) = . 



(38) 



d (2//x) 



Eq. (38) has infinitely many solutions \i = /i(A), and one needs to further 
minimize the average cost in Eq. (2) versus the average photon number N 
parameterized by A and /x = /i(A) 



In this way one can find the normalized and finite-energy states that achieve 
the minimum cost for the optimal POM. 

The solution (36) of the recursive relation (35) has some similarity with the so- 
lution for the minimum phase- uncertainty states of a single-mode field [14,15]. 
The proof of convergence of the series in Eq. (37) can be found in Ref. [15]. 
However, the matching condition (38) (instead of the vanishing condition for 
h n with n < for one mode) makes the two-mode phase estimation prob- 
lem more difficult, since one cannot exploit the properties of the zeros of the 
Bessel functions in an asymptotic approximation, as done in Ref. [16] for the 
single-mode case. 
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